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The Majumdar-Papapetrou system is the subset of the Einstein-Maxwell-charged dust matter 
theory, when the charge of each particle is equal to its mass. Solutions for this system are less 
CN ■ difficult to find, in general one does not need even to impose any spatial symmetry a priori. 

For instance, any number of extreme Reissner-Nordstrom solutions (which in vacuum reduce to 
l/^ \ extreme Reissner-Nordstrom black holes) located at will is a solution. In matter one can also find 

^O ' solutions with some ease. Here we find an exact solution of the Majumdar-Papapetrou system, 

a spherically symmetric charged thick shell, with mass m, outer radius Tq, and inner radius n. 
T^lj- I This solution consists of three regions, an inner Minkowski region, a middle region with extreme 

^^ ■ charged dust matter, and an outer Reissner-Nordstrom region. The matching of the regions, 

O , ■ obeying the usual junction conditions for boundary surfaces, is continuous. For vanishing inner 

radius, one obtains a Bonnor star, whereas for vanishing thickness, one obtains an infinitesimally 
thin shell. For sufficiently high mass of the thick shell or sufficiently small outer radius, it forms 
an extreme Reissner-Nordstrom quasi-black hole, i.e., a star whose gravitational properties are 
virtually indistinguishable from a true extreme black hole. This quasi-black hole has no hair 
^ ' and has a naked horizon, meaning that the Riemann tensor at the horizon on an infalling probe 

j^ ■ diverges. At the critical value, when the mass is equal to the outer radius, m = Tq, there is no 

smooth manifold. Above the critical value when m > r^ there is no solution, the shell collapses 
into a singularity. Systems with m < r^ are neutrally stable. Many of these properties are similar 
to those of gravitational monopoles. 



bJO. 



1. Introduction 

One can couple Newtonian gravitation to Coulomb electric fields without much effort. A particle 
with mass mi and charge gi is certainly a solution of the system. If gi = mi , then one can 
put a second particle with g2 = "^2; and gi and g2 having the same sign, anywhere. This is 
also a solution since the Coulomb repulsion compensates exactly for Newtonian attractive force. 
And of course, one can then put any number of particles, such that for each one the mass equals 
the charge and the charges have all the same sign. One can go further, distribute the particles 
continuously to make a fluid, such that the energy density p and the charge density pe of the 
charged fluid obey pe = p, and find this is going to be a solution of the Newtonian- Coulomb 
gravitation. There are other solutions for the Newtonian- Coulomb gravitation, which do not 
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have pe = p, but in these cases the density profile is singular somewhere, or one has to add some 
pressure to the fluid. 

The general relativistic version of the Newtonian- Coulomb gravitation is the Einstein-Maxwell 
system which admits a vacuum spherically symmetric, particle like solution, given by the Reissner- 
Nordstrom metric ^I21- For pure vacuum this solution represents a charged black hole when the 
mass is greater than the charge, an extreme charged black hole when the mass is equal to the 
charge, and a naked singularity when the mass is smaller than the charge. Weyl in 1917 also had 
the idea of studying vacuum general relativity and electromagnetism together jHj . Upon further 
imposing axial symmetry he showed that if there is a functional relationship between the gu metric 
component and the electric potential ip , that relation should be of the form gu = 1 + Bip + ip'^ . In 
1947 Majumdar [l] showed further that this relation holds for no spatial symmetry at all, axial 
symmetry being a particular case. Moreover, he showed that if gu = (1 + C ip)"^ , i.e., the above 
relation is a perfect square, then the spatial metric is conformal to the flat metric, the equations 
of the system reduce to a master equation of the Poisson type, and the system needs not to be 
a vacuum system, one can have charge matter with the charge density pe equal to the energy 
density p, p^ = p, thus flnding the general relativity counterparts to the Newtonian- Coulomb 
gravitational solutions. Papapetrou explored the same set of ideas, his paper was submitted one 
month before that of Majumdar Ej. In vacuum, p^ = p = 0, the Poisson type equation reduces to 
the Laplace equation, and any number of extreme Reissner-Nordstrom black holes located at will 
is a solution, as was elucidated by Hartle and Hawking [H]. For non vacuum systems, systems with 
Pe = p, the Poisson equation to be solved is more complicated and one usually imposes additional 
symmetries. Along this line, new solutions were found and discussed by Das |7] and several other 
authors P]-^]. Bonnor and collaborators, in a series of papers, have further developed the fleld, 
by constructing spherical symmetric Majumdar-Papapetrou charged stars |in]-|2SI, as well as 
stars of other shapes, ellipsoidal for instance j21]. Further developments appeared in Lemos and 
Weinberg [23] where the causal structure and physical properties of new Majumdar-Papapetrou 
stars were analyzed. All these solutions may have some astrophysical appeal, since extreme dust 
matter can be realized in nature with slightly ionized dust (for instance, a dust particle with 
10^^ neutrons and 1 proton, or an equivalent conflguration, realizes extreme dust matter). The 
equilibrium needed for the dust to be ionized, such that Pe = p always, is certainly precarious but 
possible. Note further that the Majumdar-Papapetrou solutions, both in vacuum and in matter, 
are also of great importance in supergravity and superstring theories since, when embedded in 
these extended gravities, they are supersymmetric PB|I57], and saturate a bound on the charge, 
the charge is equal to the mass, being then called BPS (Bogomolnyi-Prasad-Sommerfleld) objects. 
For instance, low-energy superstring theory admits a series of Majumdar-Papapetrou BPS type 
solutions in higher dimensions, extreme black extended objects, which can be constructed from 
extreme black holes by arraying in speciflc directions (see, e.g, ^M,), a process similar to flnding 
charged dust solutions. 

Drawing on the work of Bonnor and collaborators [201-122] and on the techniques developed in 
|25j . we report here on a new spherical symmetric Majumdar-Papapetrou star, a thick shell, made 
of three regions: an interior flat part described by the Minkowski metric up to radius n , a middle 
matter part, from r[ to Tq, described by extreme matter p = pc with an appropriate metric, and 
an exterior region, from Tq to inflnity, described by the extreme Reissner-Nordstrom solution. 
The appearance of the three regions is novel. In section II we present the Majumdar-Papapetrou 
system in isotropic coordinates, which are the ones appropriate to flnd the solution. In section 
III we present the solution in Schwarzschild coordinates to further study the properties of the 
system. In section IV we conclude. 
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2. The thick shell solution in isotropic coordinates: Equations and 
solution 

2.1. Equations and fields 

We will study Einstein-Maxwell theory coupled to charged dust (charged matter with zero pres- 
sure). The equation for the Einstein-Maxwell-charged dust system is given by (G = c = 1) 

Gai, = Sn{T^r + TT), (1) 

where Gab is the Einstein tensor. T^^^'^*' is the dust part of the stress-energy tensor 

T^r = P^aUt, (2) 

with p being the energy-density, and Ua is the four- velocity of the fluid. T^™ is the electromagnetic 
part of the stress-energy tensor 

TX = ^ {f:F,^ - IgabF'^'F^,) , (3) 

with Fab being the Maxwell tensor. The two Maxwell equations, are 

^a% = 4vrj,, (4) 

Flab;c] = . (5) 

where ja = PeUa is the electric four-current, pe is the electric charge density of the dust, and [] 
denotes anti-symmetrization. Equation (jS)) permits to define a vector potential Aa such that 

Fab = Ah,a - Aa,b ■ (6) 

Now, for a static pure electric system one can make the choice 

«a = |-, Aa = 5.V, (7) 

where U and (p are functions of the spatial coordinates, U being identified with the gravitational 
potential, and cp with the electric potential. Furthermore, it was then shown by Majumdar |3] 
in a very elegant paper, that in the special case of extreme dust matter, i.e., when the energy 
density is equal to the charge density, 

Pe = P, (8) 

the metric can be put in form 

ds^ = -^ + U' (dx' + dy^ + dz^) , (9) 

where (t, x, y, z) are called isotropic coordinates. The Einstein-Maxwell-extreme dust matter 
system ([T]) then can be reduced to the following system of equations, 

V2f/ = -47rpf/3, (10) 

¥^ = 1-^, (11) 
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with V^ 



dx'^ 



+ 



dy^ 



+ 



dz^ 



is the flat Laplacian operator. We should comment that there is 



a possible choice for the sign of the charge density, strictly speaking pe = ±p. Above we have 
chosen the plus sign, had we chose the minus sign the sign of the electric potential also changes, 
i.e., (/? = 1 + 1/U . In order to not carry a ± label we stick to the plus sign of the charge density, 
knowing the results will not be altered. 

Thus, from equation (jlOj) one sees that solutions with no spatial symmetry at all, are so- 
lutions of this system. This is because particles with charge equal to mass, or matter systems 
with pe = p, exert zero force upon each other, or within themselves, a result valid in both New- 
tonian gravitation and general relativity, so that they can be distributed at will. Solutions of 
this Einstein-Maxwell-extreme dust matter system are generically called Majumdar-Papapetrou 
solutions IHEj, a particular instance of these are the Bonnor stars |19]-|24 



2.2. Thick shell solution with three regions: interior flat, middle charged dust, 
exterior extreme Reissner-Nordstrom 

Here we want to flnd a spherical symmetric solution in which case the line element (jHl) should be 
written as 






dR^ + R^ {dO'^ + sin^ d , 



(12) 



where U = U{R), and the three-dimensional Laplacian is now V^f/ = -m-Sn (-R^fp) • The fleld 
equation (fTUj) can be solved by guessing a potential U , and then flnding p. The solution is then 
complete because pe and (p follow directly. If the solution is physically acceptable, then it is of 
interest. One such spherical symmetric solution has been found by Bonnor and Wickramasuriya 
|2IIj-|22], describing a matter conflguration from R = to R = Ro, and then vacuum outside, 
given by the extreme Reissner-Nordstrom metric. We flnd here a non trivial extension of this 
Bonnor star solution by adding a region internal to the matter that is flat, described by the 
Minkowski metric. We display for the flrst time a global Majumdar-Papapetrou exact solution 
that is composed of three regions, the interior vacuum Minkowski, the middle matter region, and 
the exterior extreme Reissner-Nordstrom electrovacuum region. 

The system is characterized by its ADM mass m, and its inner and outer isotropic radii Ri 
and -Ro, respectively. The solution for f/(r) (and thus for the system) is 






0<R<Ri 



(13) 



m 



L^ matter ' u 



1 + 



(i?o — -Ri) — (-R — -Ri) 
2-Ro (-Ro — -Ri) 



2n 



R[ ^ R ^ Rq , 



(14) 



U, 



m 



RN 



1 + ^, Ro<R. 

ri 



(15) 



With U we can flnd through Equation (fT^ the metric for the whole spacetime, and through 
Equation (fTTj) the electric potential ip. The energy density is zero outside matter, and in the 
matter it is found using Equation (jTU)), to be 



P 



3m 



Ro 



1 _ lEi 
3 R 



4:71 R^Ro 



Ri Ui 



(16) 



matter 
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with the charge density being then pc = p. Since p > throughout the spacetime, one can 
say that U given in Equations (|T!^ - (fTH|) yields so far a physically acceptable solution. Regular 
boundary conditions at the center and at infinity are obeyed, at the center the metric and matter 
fields are regular, and at infinity spacetime is asymptotically fiat, the electric field decays as 
a Coulomb field, and the extreme dust matter density is zero. Moreover, since there are two 
interfaces one has to check whether the junction conditions are obeyed at these boundaries. 
Indeed, it is found that they are obeyed, with no need for thin shells, see Section 13.2.1 for more 
details. If one puts Ri = one gets the Bonnor solution f^ with two regions, the matter region 
and the exterior extreme Reissner-Nordstrom region. Note also that on using (fTSj) on (fT^ one 
clearly gets the Minkowski metric, but not in the usual form. To do so one has to define a i 
such that i = t/t/flat? and an R such that R = f/fjat-R- The coordinates in the matter and 
electrovacuum regions would change accordingly. We will not do it here, since it is not necessary 
and gets cumbersome. 

The use of isotropic coordinates is good for finding solutions, as Majumdar has shown |1], 
but the interpretation of the solution is best done in Schwarzschild coordinates, where the radial 
coordinate r yields a measure of circumferences, or of areas, of the type we are used to in fiat 
space. 

3. The thick shell solution in Schwarzschild coordinates: Physical 
features 

3.1. Equations and fields 

In Schwarzschild coordinates we write the metric in the form 

ds'^ = -B{r) dt^ + A{r) dr^ + r^ (dO'^ + sin^ 9 dcj?) , (17) 

and the electric Maxwell field is written as 

Ao = <^(r) , Ai = Q. (18) 

One can use Equations (HI)-® to find three non-trivial equations, 

^-^ = 87rrpA, (19) 

^^^' = 87rrV+^^", (20) 



r(l-- 



rVA 



rp^ 



y/AB 



^ 



-47rpe. (21) 



Upon using the Majumdar-Papapetrou condition Pe = p one should be able to find the thick shell 
solution in Schwarzschild coordinates. Of course, it is much easier to make a simple coordinate 
transformation from the solution in isotropic spherical coordinates (t, R, 6, 0) to the solution in 
Schwarzschild coordinates {t,r,6,(f)). This we will do next. 
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3.2. The thick shell solution 

In Schwarzschild coordinates, where the the metric takes the form ()17j) . one sees that the coor- 
dinate r has some physical meaning, the area of a sphere of constant r is given by 4 vr r^ . By 
inspection of Equations (fT^ and (fTTj) one finds that the relation between r and R is 



UR. 



(22) 



Then, one gets B{r) = 1/U^. From (J22D one finds dr/[l + f ^(r)] = UdR. Thus A{r) = 
1/[1 + §j^{t")]'^ ■ Note that from (j22j) one gets i? as a function of r, in general implicitly, but in 
some cases, such as here, it can be solved explicitly. Now we have to transform coordinates in 
each region. We begin with the easiest ones, the exterior and interior vacua regions, and then do 
the matter region. The three regions are schematically represented in Figure 1. 

The exterior Reissner- Nordstrom region is trivial, and we will do it first. From Equations (|TK|) 
and ip^ one gets, for the whole Reissner-Nordstrom region 



R 



r — m . 



and in particular at the boundary, one has 
Ro = To - m . 



(23) 



(24) 




Extreme 
Reissner-Nordstrom 



Figure 1 - A schematic drawing of the three regions of the thick shell solution: the interior flat region, 
the middle charged matter region, and the exterior extreme Reissner-Nordstrom region. 

For the interior fiat region one finds from Equations (J13p and (j22|) that R\ is a function of r^ 
and To, i.e., R\ = i?i(ri,ro), and we define fi = Ri{ri,ro). Then 



/i 




--1 
m 



m . 



Then again from Equations ()13|1 and ()22j) one has for the whole fiat region 
R= 



1 + 



rp-m-fi 



(25) 



(26) 



2(ro— m) 

For the matter region, from Equations (fT^ and ()22|). one has to solve a nasty cubic equation 
for R, giving i? as a function of r. We denote this function by f{r), i.e., in the matter region, 

R = f{r) , (27) 
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and give /(r) in the Appendix A. Note that /(n) = /j. 

We can now write the metric in Schwarzschild coordinates for < r < oo. For the interior 
flat region is, 



ds' 



For the matter region is 
ds^ = -[Mf dt^ + 



< r < Ti . 



(28) 



dr2 



( 



/2(r) '"[/W-/i] 



(ro-m)^ (ro — m — /i) 

Ti < r < To . 



For the exterior Reissner-Nordstrom region one has 
ds^ = - 



m\^ , r, dr 

— ] dt^ 

r 



+ 7YzW+-'K+s-'^' 



To < r . 



(29) 



(30) 



The three regions join continuously, since the metric and the extrinsic curvature of the two 
boundaries match, as we show in Appendix B. . A quick way to see this is to note that at the 
inner boundary the gtt part obviously matches, and the Qrr part also matches since in the matter 
region one has Qrr = 1 at r = ri. In the outer boundary one has f{ro)/ro = 1 — m/ro-, so the g^ 
terms match, and as well the g^r are identical at r^- One can also check that 



df{r)/r 
dr 



so 



that the first derivatives of gu match at the inner boundary. As well, J^J''/'^ ) =2(1 — -^)^, 
so that the first derivatives of gu match at the outer boundary. The derivatives of g„ do not 
match, but that is expected for a boundary surface. So g„ is C^ at both boundaries and the other 
metric functions are C^ or higher. This means that the Israel junction conditions are satisfied, 
yielding a boundary surface (see also Appendix B). With some extra effort one can smooth out 
the grr component at the boundaries, so that grr is C^ or C^, but we will not do it here. 

The electric potential is taken from Equation (fTT|) . </? = \ — jj^ but now U is seen as a function 

of r, (alternatively, y? can be taken from Lp = \ — J\gtt\ with the \gtt\ for the three regions given 
in Equations (l ^ - lfHUj) ). Thus, we have 



V 



^ 



^ 



1- 



1 1 "^ 


(l M r-o-m-fA ' 
V^ ^ 2{ro-m) J 


< r < Ti , 


To—m 




1 




n < r < ro , 


t^ matter 


/ m\ 

1 

V r J 


m 
r 


To <r. 



(31) 

(32) 
(33) 



where in Equation (|32| . t/matter(^) is now 



a 



m 



matter 



1 + 



m 



/i)' - (fir) - /O 



2n 



mj 2(ro — m) (vo — m — /i) 

The electric potential is continuous and smooth at the boundaries. 



r[ < r < Tq . 



(34) 
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The energy density can be also calculated in these coordinates. We have 
p = 0, 0<r <n, (35) 



3m 1 ~ 



2/i 



P = 1—f \2T 7T 77 ^^^ n<r <ro, (36) 

4 7r(ro - m)^ (ro - m - /;) f/matter(^)'' 

P = , ro < r , (37) 



where Umatterir) is given in (j34|l . The energy density is a step function at Tj and Tq, but this is 
no problem, jumps in the energy density are allowed. It is easy to check that the energy density 
p is everywhere positive in the thick shell, as it should for a realistic fluid description. The charge 
density follows from pe = p. 

Again, regular boundary conditions at the center and at infinity are obeyed, at the center the 
metric and matter fields are regular, and at infinity spacetime is asymptotically fiat, the electric 
field decays as a Coulomb field, and the extreme dust matter density is zero. Moreover, as we 
have just showed the junction conditions at the outer and inner boundaries are also obeyed. 

3.3. Physical properties of the thick shell 

3.3.1. Prom low gravitation configurations to extreme quasi-black holes 

The thick shell solution is given by the equations of the last subsection (f^ - (p?7j) . The solutions 
are characterized by the mass m, the outer and the internal radius, r^ and Tj. An important 
dimensionless quantity is the ratio of the mass to outer radius, 

m , , 

a = — . 38 

To 

This quantity is a measure of how general relativistic is the system. Small a means the thick 
shell is very dispersed, Newtonian gravitation might suffice. For large a general relativistic effects 
come into play and eventually at a — > Ocrit = 1 (m — >• Tq) an event horizon is about to form. 
This is displayed in Fig. 2, where we draw the dependence of the metric, electric and density 
potentials as a function of radius, for a small and Ocrit • 

The function A signals the appearance of an event horizon, when l/A passes through a zero, 
a horizon is formed, in this case it should be an extreme horizon, since 1/A gets a double zero (it 
is not strictly a double zero, since it is not smooth at Tq, but with some care, either analytically or 
numerically one can smooth it out). The function B is the redshift function, with S = meaning 
infinite redshift. It is seen from equation ()25p that for Ocrit the parameter j\ is identically zero, 
which means 5 = in the region < r < ri . This is true also in the region Tj < r < Tq , as can be 
checked from equations (|48p and (j49p for the function /(r) in the Appendix A. Hence, at Ocrit the 
redshift is infinite at the horizon and in the whole region inside. This means that in fact a true 
black hole does not form, since inside there is no smooth manifold. One refers to the location of 
the minimum in the function l/v4, when a -^ a^rw, = 1, as the quasi- horizon, and the object is a 
quasi-black hole [23]. The product (Ai?)^/^ is an important function, it gives whether the horizon 
is naked or not, see the last two paragraphs of this section for a more specific discussion on this 
point. The electric potential (/J shows that at Ocrit, the electric potential outside the horizon is 
Coulombian, the black hole has no hair. The density field, shows a very interesting feature, when 
the star is about to form a quasi-black hole, there is still equilibrium, nothing special happens, its 
surface being a quasi-horizon with radius r^, arbitrarily closed to the extreme Reissner-Nordstrom 
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black hole radius rh. Note that the figures displayed in Figure 2 correspond to Figures 2-6 of 
PH] . The curves for A, B , {ABY''^ and show similar behaviors in both cases, since for r > Tq 
they are the same functions and for r < Tq they are power law functions. 




B 0.5 




1/2 
(AB) 0.5- 




<P 0.5 
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r 


" ~ 














Figure 2 - Behaviour of the metric, electric and matter functions in terms of the radius r . Note that r* 
is To critical. 

There are two points worth noting in the meantime: (i) It is interesting that on one hand 
these thick shell Majumdar-Papapetrou stars superseed the Buchdall bounds for perfect fiuid 
matter stars [21], where the radius of the star obeys rgtar > f^^h, with rh = 2m, on the other 
hand, when they are about to form a black hole the spacetime turns into a degenerate one, so 
that there is no smooth manifold, (ii) Although we do not discuss in detail in this paper what 
happens above criticality, i.e., for a > Ocrit (m > Tq), we note that in this case one would have a 
shell of matter inside an extreme black hole at r^^ = m. Such a solution would be everywhere free 
from curvature singularities and, following a theorem by Borde [HD], this would mean that the 
topology of spacelike slices in this black hole spacetime would change from a region where they 
are noncompact to a region where they are compact, in the interior (an example of this appears 
in the Bardeen model JSI1E21 )■ In our case this in fact does not happen, there are no solutions 
with m > To, the shell collapses into a singularity. 

We now dwell further on the nakedness of the extreme quasi-horizon. A horizon is defined as 
naked when the tidal forces suffered by an infalling particle at the horizon are infinite [32]. For 
a usual black hole spacetime, such as Schwarzschild or Reissner-Nordstrom, the Riemann tensor 
Rabcd is nonsingular at the horizon, and nothing special happens to an observer or a particle 
crossing it. But for some other black holes, such as some black holes with a dilaton field and 
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the quasi-black hole studied here, the situation can be different. Due to the acceleration of the 
observer or particle, the Riemann tensor projected in a tetrad frame freely falling with it, may be 
divergently different from the Riemann tensor in the static coordinates. For the metric ()17|) one 
has that Rf§f§ , where a hat means that the quantity is evaluated on a tetrad freely falling frame, 
is given by Rf§f§ = ~y- t" Tr ~ 7 ' ^^^ similarly for Rtiti- The important function here is the 



product (ABY^"^ . When (AB) = 1, as is the case in the Schwarzschild or Reissner-Nordstrom 
black holes, nothing special happens. But, when {ABY^'^ -^ 0, one has Ri§f0 ~ ^ and thus this 
local component of the Riemann tensor at the horizon, i.e., the local component of the Riemann 
tensor as measured by a freely falling observer passing through the horizon, diverges, and the 
horizon is naked. It is also interesting to relate this divergence to the proper time an observer or 
a particle take to make a return trip in the quasi-black hole spacetime. To see this, send a massive 
particle from a large radius r, through the thick shell, so that it turns around and comes back 
Due to the staticity and spherical symmetry of the metric one can define a conserved energy 



and a conserved angular momentum per unit mass for the particle, J = t^^, E = B{r)-^. Then 
from Equation (J17p one has 

*- ^£1!^^^. (3,) 



^2 - 5 ( 4 



dr 

We are interested in the quasi-horizon configuration, the one which is about to form a horizon but 
not quite. Call the radius of this configuration r* . Call also ^ = e where e — > at the critical 

solution. Then one has from Equation (jH^ that the particle spends a proper time inside the star 
of the order of Ar ~ ^(Ai?)^/^!^, ~ e^/^ (indeed, put r; = to simplify (this does not alter the 
result), then for 1 — — = | and — — 1 = ae, with a a number of order one and e << 1, one has 

' ' -mom' ' 

/(r) = (r„-m)(l-(e)i/2/3), and thus l/A^'^ = I - {m/r)[f{rf /{r^- mf] = e^'^ , so l/A = e, 
and B^/'^ = f{r)/r = ae, so that (ABY^"^ = ae^^^ ). For the quasi-horizon this time is extremely 
short, and at the critical configuration is zero, since (ABY^'^ vanishes for r < r^,. In addition, 
one finds Rfg^g ~ ,^^2 • Since A r goes to zero at the critical solution, the local Riemann tensor 
diverges, as previously stated. Another interesting time, is the time A t the particle takes in its 
trip for a coordinate observer, given by A t ~ -jf^ . At the quasi-critical solution this time is 
arbitrarily large. Thus, for an external observer the particle takes an arbitrarily quasi infinite 
time to come back. The observer cannot probe the inside, and associates it with an entropy for 
the star (see jHl] for a discussion on quasi-black hole entropy). 

3.3.2. Penrose diagrams 

One can draw the Penrose diagrams for the two configurations of interest, the quasi-black hole 
with a — >■ m/ro quasi-critical, and the black hole with the same mass parameter, see Figure 
3. For the quasi-black hole, the Penrose diagram is the same as for a small. In these cases, the 
spacetime is asymptotically Minkowskian, and the origin of the coordinates r = and the surface 
of the thick shell star Tq = r^, are timelike surfaces. On the other hand, for the true black hole, 
spacetime is still asymptotically Minkowskian, with rh being a true horizon, whereas r = is 
now a curvature singularity. 

Note that at criticality, {m = Vo), there is no Penrose diagram for the quasi-black hole. Since 
time has disappeared from spacetime, the solution is now degenerate. Above criticality, (m > Tq), 
there are no static solutions. 
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'■=0 




(a) quasi-black hole (b) extreme black hole 

Figure 3 - Penrose diagrams for (a) the quasi-black hole, and (b) the extreme black hole, are shown. 

3.4. Two limits: stars and thin shells 

There are two interesting limits of the Majumdar-Papapetrou thick shell solution. 



3.4.1. Bonnor stars 

If one takes the limit of a vanishing inner radius r^ = one obtains a Bonnor star [2n]-|22]- This 
solution has two regions, the interior from r = to r = Tq , made of extreme charged matter, and 
the exterior, an extreme Reissner-Nordstrom region. The Bonnor star has the same properties of 
the thick shell solution, when the star is compact enough it forms an extreme Reissner-Nordstrom 
quasi- black hole with a naked horizon (see J2S1 for more details) . 



3.4.2. Infinitesimally thin shells 

Another important limit is the infinitesimally thin shell limit, when the thickness of the shell 
goes to zero, r^ ^ Tq. The important quantity to keep track is the energy density p (an energy 
per unit volume) that must go over into a surface energy density a (an energy per unit area). 

From Equation (jHBj) one obtains, when r^ ^>- r^, p = j^ ( ^ _''" ) , in first order in Tq — r; . The 



47rr^ 



ro— n 



proper surface energy density of the thin shell should be defined as cr = limr„_>j.. ( \°_^ ) p, where 

the factor (ro — Tj) is as in Newtonian theory, and the factor 1/(1 — — ) takes care of the proper 
length (see metric (HZj)). Thus, we find 

m 



a 



Aurl 



(40) 



as it is expected for an infinitesimally thin shell. The three regions are now, an inner fiat region, 
a thin shell, and an exterior Reissner-Nordstrom vacuum. The metric is now 



ds' 



ds' 



l_!^j dt^ + dr' + r'[de' + sm^d 
dt^ + r^ (de'^ + sin^ 9 d(, 



< r < To 



1- 



(41) 
(42) 



12 A. Kleher, J. P. S. Lemos and V. T. Zanchin 

ds' = -(l~-y dt'+ '^^^ + r^ (de^ + sin^ e d<p') , r,<r. (43) 



r , 



Note that the grr component has a step at tq from 1 to 1/ f 1 — ^ j , meaning that the junction 
conditions are not satisfied there. This is a minor problem, one can always see it as a very thin 
shell (not infinitesimally thin), where there is no discontinuity, although the slope of the Qrr 



function is very high. The electromagnetic field can be taken from ip = 1 — J\gtt\ in equations 
Mip - ipSj) . The thin shell has analogous properties to the thick shell solution. When it is compact 
enough, for a -^ Ocrit = 1 (i-e., rra — >■ Tq) it forms an extreme Reissner-Nordstrom quasi-black 
hole with a naked horizon. When the precise equality holds, a = Ocrit, the redshift function B 
is zero not only at the horizon but also in the whole region inside, meaning that in fact a true 
black hole does not form, since inside there is no smooth manifold. This signals either a change in 
topology or instability. Indeed, one can ask what happens for a > Ocrit , "m > Tq ? In such a case 
a solution does not exist, the thin shell collapses into a singularity. Note also that for the same 
mass parameter m there is also the branch of pure extreme Reissner-Nordstrom black holes. 

4. Conclusions 

We have found an exact solution of the Majumdar-Papapetrou system, the thick shell star so- 
lution, consisting of three regions, an inner Minkowski region, a middle region with extreme 
charged matter, and an outer Reissner-Nordstrom region. The system is neutrally stable, as all 
the systems of Majumdar-Papapetrou type are. For sufficiently high mass, or sufficiently small 
outer radius, at almost the critical value, the thick shell forms an extreme Reissner-Nordstrom 
quasi-black hole, with no hair and with a naked horizon, i.e., the Riemann tensor at the horizon 
on an infalling probe diverges. At the critical value there is no smooth manifold. Above the 
critical value when m > Tq one has an extreme shell inside an event horizon collapsing into a 
singularity. All these properties are similar to the properties found for gravitational monopoles 
when gravity takes care of the system [HI]. In another work we will explore these similarities 
which are much too striking to pass without some attention fl^ IHH] . 

Thick shell solutions are rare in the literature. For a discussion of the generalization of the 
Israel thin shell junction conditions [SHI to thick shells see ^37j, and for an application to the 
covariant entropy bound in Tolman-Bondi spacetimes see pS] . 
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A The isotropic coordinate R as a function of the Schwarzschild 
coordinate r in the matter region 



We want to invert the equation 

r = -Rt^matter = R 



Ro \ 2Ro{Ro — Ri) 



(44) 



In order to condense expressions, define 

E 1 



f-3«- '=f + 



2E 2Rf\ „ 
+ ^] Ri 



3F 27 



(45) 



with 



Ro \ 2i?o 



m 



(46) 



2Rl{Ro - Ri) ' 
where Ro = r^ — m as in (j24j) . and R^ = Ri{ri, Tq) as in (j25p . Then through further definitions 



A 



b b^ a3 



\-2 



Ib^ a? 

~A^2^' 



(47) 



the solution is i? = /(r), with 

/(r) = ^/.-^^— ^-1— ^V-3. 



2 A+B A-B 

3 2 2 



For r and r\ small one finds 

fi ^ r 2(ro - m) 
2ro + m 

where we have used (j^ and (P3jl . 



(4J 



(49) 



B The matching conditions of the the thick shell 

We will follow Israel [Sn] to show that the two boundaries at Tq and r^ are boundary surfaces, 
i.e., do not need extraneous thin shells. We need to show that the metric, or first fundamental 
form, is continuous at the surface, (7^5 = (?~^, and the extrinsic curvature, or second fundamental 
form, is also continuous, K^^^ = K~f^. For the extrinsic curvature, we note first that, since the 
surface we use is spherical the normal Ua to the surface has only a radial component Ur = ^Jg7r- 
Then the extrinsic curvature has the form 



Kl 



-'^r^ afi 



c?x" dx^ 



di^ di^ 



(50) 



where ^* are the intrinsic coordinates of the surface. Now we analyze the junction at the outer 
and inner surfaces. 
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Outer boundary, r = Tq : For the surface Vo we adopt the metric 



ds^ = -il- — \ dr^ + r\de'^ + sin^ d(l)'^) , (51) 

with the intrinsic coordinates ^* being ^* = {t,6,(J)). Then one has /(ro)/ro = 1 —m/ro, 
so from Equations (J29|l and (J30|l one sees that the g^ and g^ match. The g^^ and g~j. also 
match, as well as the terms for the angular part. Now, at Tq coming for the interior one finds 
K:;^ = -rirVit ^1;. By construction ^ = 1- One also has n^ = ^ ,-2 ,i(/_/,) U = 

r (ro-m)2(ro-m-/i) 

^. Noting that f |.„ = 1 one finds r^ = '"^^°7'"^' • Thus, iT" = !Zi(l2^. The extrinsic 

ro 00 

curvature coming form the exterior Reissner-Nordstrom region can also be calculated, one 
finds i^^ = "^^^°7 • So they match. One can also check that K00 and K^pj, also match. 

' o 

Inner boundary, r = n : For the surface ri we adopt the metric 

ds^ = -(iL\ dr'' + r\de^ + sin^d(f'). (52) 

From Equations (|28p and (j29p one sees that the gtt part matches, and the grr part also 
matches since in the matter region one has ^f^r = 1 at r = r;. Likewise for ggg and 



Now, at Ti coming for the matter one finds K^^ = —UrT^tt ^^- By construction ^ = 1. 
One also has n^ = 1. Noting that ^1^; = ^ one finds F''^ = 0. Thus, K^^ = 0. The 
extrinsic curvature coming form the interior fiat region can also be calculated, one finds 
K~^ = 0. So they match. One can also check that Kgo and K^^f, also match. 

So the metric and the extrinsic curvature are continuous as is required. The electric potential ip 
and its first derivative at the boundaries are continuous as is also required. 
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